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VI. On the pear-shaped Figure of Equilibrium of a Rotating Mass of 
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Introduction. 

This is the sequel to a previous paper on " Ellipsoidal Harmonic Analysis " (' Phil. 
Trans../ A, vol. 197, pp. 461-557). I here make use of the methods of that paper, 
and for brevity shall refer to it as " Harmonics." 

The sections 1 to 4 are preparatory, and might have been included in " Harmonics," 
but seem, more appropriate here. Section 5 is an independent investigation of so 
much of M. Poincake's celebrated memoir on rotating liquid* as relates to the 
immediate object in view. 

It is not necessary to say more here, since I give a short summary in the last 
section. 

§ 1. The Harmonics of the Third Degree. 

It was remarked in § 7 of " Harmonics " that this group of harmonics may be 
determined in rigorous algebraic form ; I now proceed to do this. 

The requisite formulae are contained in § 6 of " Harmonics," and I take the several 
harmonics in succession. Throughout this section we have, of course, i = 3. 

,5 = 0; type OEC, and %(v) = q P,(v) + fa % P£{v\ with q = 1. 
The equation for j3a is 

P(T — 4 + #r 1 + i/3<r ' 

If we write 

(5 1 )»=l+15/8 B , 

the proper solution of the quadratic equation is 

/3cr = 2{B X - 1). 

* < Acta Mathematical vol. 7, 1885. 
(305.) 17.3.1902. 
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Also 



Then since 
we find 



a \ 1 

P 3 (z>) = f ^ 3 — f v, P P / 2 = 15v (?^ 3 — 1), 

f M - jjjg ( A - 1 + 5£) v (v s — ^— J- j . .... (1). 



* = 1 ; type OOC, and P*(v) = \/^ -'? ^'^ (") + #&W („) ), with g/ = 1. 
The equation for /3<r is 

/3<r + P . 3 . 4 = ^S.i^ 3 > , 

or /3<x + 6/3 = ^~ ± ~Yo~ s 

■J ~t™ a /3CF 

If we write 

the proper solution of the quadratic equation is 

(So- = 4 (B, 2 — 1 — J-/3), 

XX-J-OO 

,. _ _-" w-«~» ___________ ■■ 1 _„. . Q TY fl ~ZZ—~ txnm** -f _ TT r 1 I h /"¥ «•»»»»■■» I 

ffl ~ 8 + /So- ~ 4(i? 2 + 1 - W dna ?1 ' ™ d £' Wltft & " X> 
Then since Pg 1 (v) = -f(5v 3 — 1) (^ 3 — 1)*, P 3 3 (v) = 15 (v 9> — 1) (v 2 — l) 4 , we find 



5=1, type OOS and ^g 1 (z/) = g^ 1 (v) + fSqzP* {v) 9 with g^ = 1. 
The equation for /3cr is 



l$0* —mm 2"yQ ,0,4: 



■^/j^-^O, Zjj "JOj Oj 



2/ ^* * 8 + #0- 

A comparison with the last case shows that we have only to change the sign of /3 ; 
accordingly if 

(P 8 ) 3 =l+t£(l + y3), 



we have /3<x = 4 (J5 3 — 1 + f/3) 

^<9q X 



— — — —j-^ — -— -•-- — "—— with 0*1 •■— 1 
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On substitution we find 



£/'•**>. - ■_ .~-\. / o .| \i / o O "j" fj 4x>g 



^( v )^~(B s -l + 3p)(^-lf( v 



5(1-/3) 



• ■ • (3). 



s = 2 ; type OEC, $ 3 3 (v) 
The equation for /3cr is 

/3cr 



% P 3 (") + % F z (")» with q % — I 



M!f 3 ' !} i 3 ' 2 } - 15/3 3 



4 — /3<x 



1 — l/So- 



We have already defined {B^f as 1 + 15/3 3 , and find the proper solution of the 
quadratic equation to be 

f$ar = — 2 (B 1 — l). 



A ISO 



2q 



o 



IMLiMi 

4 — /Scr 



4(ff x - 1) 



, with q % = 1. 



With the known values of i\ and of PA we find 



8 



3 > 



* a '< v > = T (1 + 3/S 



i? 1 ) z/ { ir — ... 



« • • . « I te #» 



5(1-/3), 



A comparison with (l) for s = shows that the last factors in each only differ in 
the sign of B v 



/ 3 _ ^ 



5=2; type OES, P 3 5 v - ; — . v 3 _ 1 
Since P 3 3 (x') =15^ (y 3 — 1), we have at once 



Pi{v) = 15^ (^- 1)* v 



,a 






» • » 



• • . I }* 



$ 



3 ; type OOC, P 8 3 („) = a/^~JeI \fiq^{v) + <Z S 'P 3 3 ("B with ?/ = L 

V 3 -j 



The equation for /3<r is 

y3tr = 
We have already denned 



|/3 3 {3, 2}{3,3} 

8 — /3er — 6/3 



- 15^_ 
8 -- /3<x - 6/3 



TO 2 - 1 - 1/8 (1 - /3) 



and find for the proper solution 



/3<r=4(l-P~ P 3 ). 
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Also 



2ft 

~/3 



a 



nd 



- {3 ? 2} {3, 3} 

8 - /3<r - 6^ 



6 



/3 £ 



(xA> '"""'- 1 ~j-" xP) 5 



ft 



ft 



-. = ■ "-, with qJ = 1. 
ft 3(/„ 



Whence on substitution 



20 



P» 8 (") = f(l -**)("» -i 



l. + # 



ft 



v 



2 



3 - /3 + 2B S \ 
5 (1 - /3) ",, 



9 * ♦ 



(6). 



,s' = 3 ; OOS, :p 3 3 (i/) = /3^ x P3 1 (p) + ^3 A 3 ( v )> "v^ith. § 8 = 1. 
The equation for /3cr is 



f3or = 



j/3 »{3, 2} {3, 3} 
8 - /So- + 6/3 



We may derive the result from the last case by introducing 

W=l +1/8(1 +0), 
and changing the sign of /J, so that 



/So- = 4(1 +1/8-^ 



2 ?i 
Whence on substitution 



16 



/S 



(5 8 — 1 - f£) } w ith q a = 1. 



WW — (l+P-BJV- If (f - ~^z^] 



» s 



(7). 



The forms of the corresponding functions of fi are the same, except that (1 — /x 2 )^ 



and 



1 + /3 
U-/8 



/rM replace the corresponding factors, 



I have not determined the cosine- and sine-functions, because they may be written 
down at once from the results already obtained. The three roots of the fundamental 

cubic are j> 3 , /x s , and — ~— . Hence we have only to replace v % by this last 

function in the seven formulae (1)-— (7) in order to obtain functions proportional to 
the seven cosine- and sine-functions. If the definition of the latter functions is to 
agree with that given in " Harmonics/' the factors must be determined appropriately, 
but the question as to the value of the factor will not arise here, 
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§ 2. Change of Notation. 

It will be convenient, with a view to future work, to change the notation, and I 
desire to adopt a notation which shall not only agree in the main with that used 
in "Harmonics," but shall also facilitate reference to a previous paper on Jacobi's 
ellipsoid (' Roy. Soc. Proc., ? vol. 41, pp. 319-336). 

I write 



K 



1 +"/3 



K" = 1 



K 



It may be noted that what I here write k was denoted by k' in " Harmonics," and 
vice versa. 

I have in general written the current co-ordinates v 9 fi, c/>, and the ellipsoid of 
reference v , so that the squares of the semi-axes are 



W f v 2 



1 + ff 

1 - fi) ' 



*»(„ »- i), *v 



I now propose to write as the squares of three semi-axes of the ellipsoid of 

reference 

e 3 cos 3 y, c 3 (1 — /c 3 sin 3 y), e 3 . 

Comparing these two we see that 

h = ck sin y, and v Q = 



tc sm 7 



For the current co-ordinates I retain </> and write 



z> 



# sili ^ 



? jjl = sin- ft 



The three roots of the fundamental cubic are therefore 



v 



2 



AT 8111 " ^ 



/X 



Sir 0, ~ — - T - = - (1 — fc " COS" r/>). 

1 — /} « 3 v ; ; 



Hie rectangular co-ordinates x 9 y, z are therefore now expressible as follows :-— 



X = 


c sin. 7 

sill -*|r 


?/ = 


c sin 7 
sin -^r 




c sin 7 


2 = 


~» 



VOL. CXCVIII. — A, 



sm -^ 



cos i// (1 — k 3 sin 3 #)~ cos <j> 9 



. (1 — k 3 sin 3 \jj) h cos # sin <£, 



- . sin 9 ( 1 — k /3 cos 3 <£)* 






(oj. 



/ 



2 R 
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These give 



(K? 



+ 



y 



cos 3 ^ 1 — « 2 sin 3 i|r 



+ z 



9, 



• o 

c, u sin" 7 
sin 3 i|r 



At the surface i/r = y, and we have 



'7 1 Qi 

cos^ 7 1 — » /r sin" 7 



3 _ c % 



In the formulae for the third harmonics, in every case but one, and in two out of 
the five harmonics of the second degree, there occurs a factor of the form {y 1 — con- 
stant) ; in each such case I write that constant in the form 5 3 /k 3 , and q f ° 2 = 1 — (f. 
Thus q will have a different value for each harmonic, 

It has been already remarked that for most purposes it is immaterial by what 
constants the several functions are multiplied. Although it would be easy to 
determine the constant in each case so as to make the function agree with its value as 
defined in " Harmonics," yet I shall not take that course, and shall omit factors as 
being in most cases redundant. 

For the sake of completeness I will give the first and second harmonics in the new 
notation, as well as the third. 

Since the harmonics of the first degree are expressed by 



f i (v) 



v, 



Pi 1 (y) 



V 



2 



\ 



1 + P$ 



l/W = (^i)= 



it is clear that in the new notation 



f i (/*) 

o 1 (*) 



sin ^ 
: sin 6 } 

= (1 — K' Z COS 2 </>)% 



Pj 1 (v) = cot \fj, 



Pa 1 fa) 

Gl 1 (*) 



w w 



(1 - * 3 sin 3 ^)» 



sin ^ 



(1 -K*sin 3 0)*, f^/x) 
cos <j&, S] 1 (<£) 



cos 

sin r/> 



(9) 



It appears from § 12 of " Harmonics" that 



f 3 (v) = V 3 + 



a 



IPs 3 ( V ) = ^ H~ 7 



where 



7 _ 7? -. 2 7' 






and .B 3 = 1 ~f 



In accordance with the notation suggested above, let 



r 



2 + ^ 



« 2 "~" 3 (1 - /9) 
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Then substituting ; -; for /?, we ii 



and for both cases 



Hence 



9: 



3 



1 + /c- =F (1 — kV 2 )*} 



K 



2 



2 



2 



2 — Sg 8 
1 — 22 



<>/>2 ' 



U 3 (") and ^(„) 
$ 3 fc) and ^ (/t) 



1 — g 3 si 11 '' ^ 
sin 3 i^ 



/c 



1 — - sin 2 8 



9: 



tf' 3 



S 3 (<l>) and at/ (<£)=!- -y 3 cos 3 <£ 



• » a » * 



(io) : 



2 3^3 

are k 2 = cf t—~9% > an< ^ 2 s = i [1 + k" ^ (1 — kV 2 )*], with upper sign for the 

first and the lower sign for the second. * 



It appears from (19) and (20), § 7, of " Harmonics" that 



P. 1 M = 

p* 1 0*) = 

C, 1 (*) = 



COS^r 

sin 3 ^r J 
sin 0(1 
cos <£(1 



k % sin 2 #)% 

K 2 COS 2 (j>) h 



and from (21) and (22) that 



s, 1 w 



(1 — # 3 sin 3 i|r)* 
snr ifr 

= sin cos #, 

= sin <£ (1 — K r2 cos 2 <£)* 



Lastly, from (25) and (26) 



p» a w 



COS \/r (1 — # 3 sill 3 ^r)* 



Sill" y 

p/( /A ) = cos6'(l -£&ux*0)\ 
jS^ 3 (</>) = sin <j> cos <£ 



« 



• (11): 



I • » e 



• ( 12 )- 



» « • 9 



(13). 



Turning to the harmonics of the third degree, we found that in the' two cases where 

the type is OEC, 

4 + B l \ 



i, („) and $ 3 2 (v) = v\r 



5(1 -0), 



«U XV -u 
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If we put 



q* 4 + ,/ij 



AT 



5 (1 — /8) 



, we find 



f/ = f [1 + K S T (1 — I -K 2 + ^ 4 ) } ] 



and 






AC' 



o. 



T 



4 — 5<2 
* 3^4^ 



Therefore, with the above alternative form for q 



$, („) and f 3 2 (v) 



1 — <f sin 3 i/r 
sin 3 1^ 



f 3 (/.) and $3* (/*) 
C 3 (</>) and 0/ ((/>) 



sin # ( 1 



K e 



( f 



m\ % 6 



fC 



j* 



I •—— cos 3 (/>)(! — k 2 cos 3 <£)* 



, (14). 



Again in the two cases where the type is OOC we found 



P, 1 (v) and P 3 3 (v) = v 



.t o — /3 + ^K$\ 



5(1-/3) A 



P 



1 4- /3\ 



Putting 



we 



find 



K 2 5 (1 — >S) ? 

<2 3 = | (1 + 2/c ^ (1 — /c 3 + 4/c 4 )*) 



ct/XXvA 



fw2 



2 9 2 ~~~ 5 r 

- y 1 - 4? 3 



Therefore, with the ahove alternative form for <f, 



P 3 * (,) and P 3 S 0*) 
P,i ( M ) and P 3 3 fo) 



cos^(l — g 3 sin 3 ^r) 
sin 3 ^ 

— (1 — k 2 sin s 0) 1 ' ( 1 — -5 sin 2 ) , 



/c 



7-x 



C3 1 (<£) and (£ 3 S ($) = cos <£ ( 1 — ' v /3 cos 2 <£ 



, 1 I u|« 



/ 



In the two cases where the type is OOS we found 



f 3 ^)andf 3 » 



9 3 + /3 + 2Z? 3 \ , u. 



5 (1 - /8) 
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Putting 



we find 



and 



K 3 ~ ""5(1 - /3)" " ' 

t/ = J- (2 + k* qp (4 - kV 2 )*) 



K" = g 






Therefore, with the above alternative form for g ! 



W W and i 3 3 (*) 

is 1 (/*) and f 3 3 0*) 



(1 — k% sin? yjr) h (1 — </ 3 sill 3 yfr) 

sin 8 i/r 



cos 0(1 — -r sin 2 ) , 
\ <2 2 



.'2 



JSg 1 (<£) and «S 3 3 (<£) = sin <j) 1 — -—■ cos 3 <£ 



# 



\ 



» » I X V* /• 



y 



The seventh of these harmonics, which is of type OES, stands by itself. We 
had 



P 3 3 (v) = v (v 2 — 1 )* I v 2 — ; 



1 + /3 V 



This gives in the new notation 



l-#,/ 



s» 8 W 



ec>ST/r(l — /c 3 sin 3 -v/r)* ? 
sin 3 -v/r 

sin cos 0(1 — • k 2 sin 2 #)* , 

sin (f> cos <£ ( 1 — k' 2 cos 2 <£)* 



» I JL $ 1 * 



The formulae (9) to (17) give the fifteen sets of three functions constituting the 
fifteen harmonic functions of the first three degrees. It would be easy, although 
somewhat tedious, to find the coefficient by which each function is to be multiplied so 
that its definition may agree with that of the previous paper. 



§ 3. Expressions for the Solid Harmonics in Rectangular Co-ordinates, 



The three roots of the original cubic equation were v 2 , /x 2 , : — „ — , and in the 
new notation the three roots of 



:tr 



f Z" ., . ., 1 a 1 — K 3 COS" (f> 

ro 3 — l//c 3 W — 1 cd- ' /r Sill" yjr K" 



3 1 
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Hence it follows that we have the identity 



'V 



4- - 

O 1 / o 1 o. 

ft)" ™ I//C" ft)" 



t 






Putting: o> 2 = 



I ftr 






c v k- sin* y = c~/r sm v y M m i:±. 11 _„/ _ js:_.„. „_-, 

' / (l/* 3 -©^! -ft) 3 )© 3 



«* 



o ? 



,i/3 



P tj > l u ii um „... ... „„„._ . 

""fj \ .■) it 

q - ht - 2- 






c 3 sin 3 7 



/:> 



2 



, + e 2 sin 2, y = r-« — ^r~TT (1 — <7 S sin 2 *//)(! — — sin 3 0) (1 — ■ —, cos 2 <£). 



This expression, together with those for x, y> z in (8), enables us to write down the 
results at once. As before, I drop the several factors as being redundant for most 
purposes. 

From (9) 






»i 1 Wfi l (^)«i 1 W = y 



»i M % 00 Ox (*) = s , P^ (v) P^ (^) d 1 (*) , 

From (10) 
- \% ( V ) fj s 0*) Cj (*) and ft» (y) ftf (/x) <£/ (*) = tf V + ^ ^ - <M 



• « » 



i r ~ </ " 



(18). 



+ c 3 g s </' 3 sin 3 y . (19). 



where 



Cf 



1 + * 3 "F (1 — kV 3 )*] , and k 3 = q* 



1 



•> > 

or/" 

4-i'f 



SO 



uoaT' 



( f ( l 
/£" — #" 



X ^*/ * 



From (11), (12), and (13) 

p 3 i (v) p 3 i m c 2 j ($) = ** , & 1 (>w o*) s, 1 w - f , p 2 - w p» 5 (/*) a* 8 (*) = *y 

(20), 



* ft 1 



From (14) 

» 3 (") \h (f*) C 3 ft) aud f 3 :: (*) % s 0*) <V ft) 



.1 . /.) 



q~(f 



I ((/% 2 + - T - : — :, y 
XJ - k" — q" 



Q 






+ c 2 </Y 2 si" 2 y) • (2i). 



where 



f / = I [1 + k 2 ^ (1 - | k 3 + k 4 ') j ] , and K a = </ : 



4^ 






4r/ 



o ? 



so 



t/iiaij 



r2 . 



K" 



( l 



2 



3 — 4q* 



From (15) 
Pa 1 W P 3 X (/*) <^ 1 (^ ^d P ;i 3 (,) P/ (/x) (^» («/,) 



o , /o 



r^/ 



x (<fx 2 + ~^^- A 2/ A 



/n. o 



i/ ^' 



+ c J (fq' 1 sin 2 7) . (22), 
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J 



where ^ = |(1 + 2k 3 =F(1-k 8 + 4k*)*) , and /c 3 = g 3 ± — %L 

ft ft 

so that 22 == 1 — 4^ 2 , 



^ 3 — # 2 



From (16) 



Pa 1 W fs 1 0*) * 3 ] (*) and fc» („) f a » (^) ^ ft) = y (g V + ^ f - qV 

K — q 



where g 2 = 1 (2 + * 3 -F (4 — kV 2 ) 1 ) and k 2 



+ c*<fq f * sin 2 y) . (23), 



1 - 2^ 3 J 



,2^/2 



so that ' o i = ~|- (I — 2g 2 ) . 



Lastly, from (17), 



P./ (*) P 3 a (m) Sa 2 '(*) = : ^ • ....... (24). 



It is easy to verify that each of these expressions satisfies Laplace's equation. 



§ 4. The Expression for the Q-functions in Elliptic Integrals. 

In this paper I drop the factors Q& and E which were found to be necessary when 
the Q-functions were expressed in series. 
We make the following definition :— 



%' ("0) ©/(*o) - [f / (*o)]* J ^ r^ (z ^ (j , 3 _ 1}J (l/3 ^±|)l 



and a similar formula holds for P/ Qj. 

It is clear that |J/ may be multiplied by any constant factor without changing 
the result ; hence we may use the forms which have been found in §§ 2, 3. 

The notation must now be changed. 

1 1 

We have v = — : — r and ^ = — r — . Therefore, when ii is adopted as variable, 

« sin y u k sin 7 r J - ? 

the limits are y to 0, and the sign of the whole is changed. 

But 1 COS yfr 7 

civ = — — .—■„-•■-- aw, 

tc Sir yjr 

and / , v / 2 x + #Y cos^(l - * 3 sin 3 ^ 

v 7 \ 1 — p/ /e~ sin" y 

Therefore p _ ^g [ y <#" 
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In accordance with the usage in elliptic integrals, I write 



A 



2 ■= 1 — . k s sin 3 



under the integral sign, or ] — k 2 sin 2 y outside the integral. 
I shall also for brevity write 



A]" = 1 — <jr" sur^i 

under the integral, or 1 — q 2 sin 2 y outside the Integral 
We have then 



m (v ) m? w = k i%- (v )j 



rv 



d^jr 



o m (vw a * 



I apply this formula successively to the several functions, as given in (9) to (17), 
and introduce the abridged notation just defined, but I do not reiterate the special 
meanings to be attached to the symbol q in each case. 

Since ^%(v) = 1, we have (dropping the now unnecessary suffix 0), 



Wi (v) €i (*) 



A 

tc p sin 3 -fr dyjr 
sin 3 7 ] A 



Pi 1 (*0 Qi 1 (*) = k cot 2 y 



7 







tan 3 yjr 

A 



a\fj, 



^ {v) ai i („) = ^. P s ^± #! 

v 7 v 7 sin" 7 J A 3 T 

f ,(v) %, (v) and » 8 »(v) «tf (v) - *£- P j^r/^ 

sm 7 j q A-j A 



P, 1 W Q, 1 (v) 



# cos 3 7 
Bin 1 7 



■v sin 1 ^ 
o cos 3 a/t A 



tty, 



W (*>) ©2 1 (l>) = -7-7- -t/ #« 



(Lj). 



i V ;V2 V"7 Bin .; J ooh 2 ^A 3 V 



ffc,(v) ® ;! (,)andf/(,) ©,*(„) = ^- f^ety, 



t^A* 

sin fi 7 J A/A 
#cos s 



P^W QaM") and P 3 3 (») Q.^) = -^- P^^ #, 
^ ' o v / -cu \ / sin r )fy fc cos 3 ^rA/A r 

Ss 1 W *V (*) and fc»<„) «,»(„) =f^ £ ~t #> 



sin 7 J cos 3 *& A" 



o rf^ 
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All these integrals are expressible in terms of the elliptic integrals 



F=[' 



7?> M>^> n= f f 1 4- 

It will, however, be found that in fact the coefficient of II vanishes in every case. 
The cases of i = and i = 1 are very simple, and we have 



$ ®o = «# 



$i©] 



fC 



Sill" 7 \ K 



2 






/ j 1 x 

P^Qi 1 = k cot 2 y ( — A tan y -— — j£ 



/e 



'3 



1/1T* 1 



W<&1 



/cA 3 
sin 2 7 



1 1 

- — /' —4 ■ I 

q -*- I o fi) a 



sin 7 cos 7\ 

/e' 3 A / " 



It is possible by direct differentiation to verify the following results, although the 
verification will be found pretty tedious. 



sin 4 ^!r t , (2 
dxfj = — 



ia«l imk 



; W _. «8 ( i _ 2(f) 2£ - 1 



22 W _ J*; 



2 2 Y 



1 ^ A sin i|r COS i/r 

27YVZ 2 1) ^ + "2YV^7)A? ' 



sin 4 ifr 7 t 

COS^ -v/rA 



^ F + "~i~7r" ^ + "^ A tan r// ? 



/r /e 2 /c' 3 "~ ! /c' 3 



f sin 4 -v|r 7 . 



2 jp 1 + /c' 3 ^ sin i/r cos ^ 



/e 



a;V 3 



*V 3 A 



sin 4 ^ 1 1 + /c 3 tan\lr ^ ox . „ 

cos 3 \^A 8 r /eV 3 A 4 # 4 A L v. « / r j 

These are all the integrals needed for the harmonics of the second degree, 
case of the first we have 



K 



2 



1 - 2? 3 * 



In the 



Thus the coefficient of II vanishes and the results are 



— *A 



P 3 (,)® 2 (^)andi/(z,)(a 3 3 (^=; 



1 — 2? 2 „ 1 — 2<f „ , (1 — 2 2 2 ) A sin 7 cos 7* 
sur 7 |_ 2r/c/°~ 2qY l "* 



2<fq'* A, 



P* 1 W Q* 1 W 



2 



# cos- 4 7 
sin 4 7 



1 7TT . 1 

..a ^ n 



fC 






A + — A tan y 



/e 



f s 1 (") € 2 X l>0 = 
P 3 3 (") Q,» (v) = 



KLX 



sm 4 7 



'2 



2 ~ 7 1 + #'* 77 sm 7 cos 7 



K 



icW 



fC fC imk 



fC cos 3 7A 
sin 4 7 



1 



2 ,,'3 



K«K 



J? \ + ^ V _L tan ^'^ — (1 + tf 2 ) sin 3 7)* 
"" ""/S? 4 " + ^A ^ 



In the first of these (f = £[1 + /c 2 =F (1 - kV 2 )*] and ac 2 = q*\ 



2 2 ~ 3 ? 



2 



1-2? 



o „ q • 
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The following integrals may also be verified by differentiation : 



f 



sin 6 -^ 



d\Jj == 



? 3 (4 - 5tf) - /e 3 (3 — 4g») 2fg'* + * 3 (3 - 4g 3 ) 

H + 



2gfi 2 , 3 («2 - 2 2) 



2/cYq'* 



F 



k 8 (3 - 2j 3 ) - 2cf 



+ f 



2^V 3 (* 3 - f) 



9 . /9 



«<n 



A sill l/r COS -^r 



..& 



(26) : 



sin 8 'xjr 



COS 3 ^ A/A 



city 



*8 (1 - 4g 3 ) - g 3 (2 



2#y 4 (a: 3 - (f) 



5q*) 4^ 2 — 

* 9^V 4 









A tan ^r 



1 +_2£*j- t£ (2 g 3 4- 1 ) 

A sin i|r COS yjr 



(27), 



sin 6 ilr 7 
7 # 



1A8 



A/A 



« 2 (1 - 2 g 2 )-g 2 (4- 



O ..9, J9, 



9 /9. 



% 3 ) „ , 2?V 3 - ** (1 - 2r/ 3 ) „ 2j¥ s +_*V S 

I O 4 /9 / 9 c >\ ■*- I o 9 '9 /•> / 9 9\9 

2tc"q q (tc — T) 2k 1 k z q z q* (/r — q 1 )-' 



Sin -\|r COS i|r 



A sin -v/r cos t|r 



sin 6 -^ 7 , 



COS 3 i/r A* 





*'»(**- ff 3)»A 22 ,a (/e a 


9 — . *- 3 

/J — — ft, y— 
— P _ 

a?V 3 


2(1 — /cV 3 ) /;7 f sin -\|r cos yjr | A tan ^ 
/r/e h fc"/c l A /c l 



v2\2A 2 



5' 



S 



(28), 



« 9 e a \ .W t/ /a 



Now in (26) we have to put 



in 



(27) 



and in (28) 



4 _ 5^2 
kt = q z 



K 



— q 



K 



T 



3 - 


- 4J 8 J 


2 - 


,-■ 9 

- sq* 


1 - 


- 42 s ' 


4 - 


- 52 3 



2<Z 3 



Introducing these values, and taking the integrals between the limits y and 0, Ave 
find : 



$ 3 <& 3 and f 3 W 



*A, 4 
sin 6 7 



2a^2 4 2 



1 



2«4/v'3 



F- 



f /9 



2q± + 5?' 



1 



' 7 / "™|— 



(4^2 — 1) A sin 7 cos 7 



22* 2 '* 



i«Jk"| 



P 3 iQ s i and P 3 3 Q 3 



/c cos 3 7A X 4 f 4^ 3 — 1 r 1 — 5r/ : 



sm u 7 



2q 



A 



9,*!Wa'± 



E 



2q i q /4> ' wq^q* 

1 — 7 g a _ (1 _ 5 g 2 .- 2g 4 ) sin 3 ry \ Atan 7 1 

U\Y % ) V" 



(30). 



. (31). 



P 3 ^3 1 and ^ 3 8®3 J 



/cA^ 3 
sin 6 7 



(1- 



2g«) ( 2 == 3g 3 ) 2 - 1 W 3 

6/c 3 2Y 4 6/c 3 /c /3 ^ /4 ' " 

5^ 2 + 6^ 4 - # 3 (2 — 11$Y 3 ) sin 3 7\ sin 7 cos 7 1 



( 0Z1 1. 
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_ kcos 3 7 A 2 f 1 + ^ 2(1 - *V 2 ) 

3 ^ 3 ~ sin 6 j { ' *V 2 ~ ' «V 4 



i A + * 2 - (1 + « -t ) sin 2 7 \ tan 7 ] 



In (30) j8 = f[l + K» T (1-1^ + «*)*], * = tf±=M. 
In (31) g» = |[1 + 2k* =F (1 - *»' + 4**)*], «» = g* ^g . 
In (32) <f = |[2 + k* T (4 - kV 2 )*], k 3 = r/ 3 ~-"|^ • 



§ 5. Bifurcation of Jacobi's Ellipsoid, 

If a mass of liquid be rotating like a rigid body about an axis, x, with uniform 
angular velocity o>, the determination of the figure of equilibrium may be treated as a 
statical problem, if the mass be subjected to a potential -|co 2 (^ 3 + £ 2 )- 

The energy lost in the concentration of a body from a condition of infinite disper- 
sion is equal to the potential of the body in its final configuration at the position of 
each molecule, multiplied by the mass of the molecule and summed throughout 
the body. In the proposed system, as rendered a statical one, it is necessary to add 
•Jco 2 (y 2 + ?}) to the gravitation potential before making the summation. If A denotes 
the moment of inertia of the body about x 9 this latter portion of the sum is ^Ac* 2 , and 
is therefore the kinetic energy of the system. 

If dm l9 dm % denote any pair of molecules and D l9j the distance between them, and 
E the energy lost, we have 

If the system had been considered as a dynamical one, the expression for the 
energy of the system, say U, would have resembled that for E y but the former of 
these terms would have presented itself with a negative sign. 

It is clear that the variation of \Ao? 9 when the moment of momentum is kept 
constant, is equal and opposite to the variation of the same function when the 
angular velocity is kept constant. 

The condition for a figure of equilibrium is that U shall be stationary for constant 
moment of momentum, or E stationary for constant co, in both cases subject to the 
condition of constancy of volume. The variations in question lead to identical results, 
and I shall proceed from the variation of E. 

Zi fc> A 
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\p 



=f 



00 



du 



(u + a 3 ) ¥ (w -f- 6 3 ) 8 (u + c 3 )* ? 



the internal potential of an ellipsoid of mass M and semi-axes a, &, c is 



4 .M 



a da b clb c do 






Hence 



X 

2 



f dm^dm^ 






00 



JO 



\Ia ULiLi l wrm 






C?77l. 



Now if A y B y C denote the principal moments of inertia of the ellipsoid about 

Xy 1^^ Zy 

| xhlm = %(C + B-A) = \Ma\ 

and similar formulae hold for the two other axes. 
Therefore 



X 

2 



dm, dm 9 

_ 



W 



x "T" T ^ "1 "4"" ~r -4- C '~;""' 



But since A P is a homogeneous function of degree — 1 in a> 5, c, the sum of the 
three differential terms is equal to — W. Hence this expression is equal to ^ -M 2 ^. 
Since 



we have 



E = -&M* 



VL/ «£»&«*< 



& 3 + c 3 



3M" 



w 



If i? be varied, whilst abc and o> are constant, it is stationary if 



da V db oM 



56 + 1 C fc + 3Jf'" 



or } 8c 

/ 



0. 



S& SZ> Sc 
— + "-- + 

abc 



0, 



Eliminating Sa ? 86, 8c we have the well-known conditions for Jacobi's ellipsoid 



9.7 9. 



2fc> 3 6_ 
3Jif 



2o> 2 c g 
~3Jf 



-1 / ^ _ ^ 

fr 3 \ ' ^a- &■ 



^ 7 <Z^ 

a da ~ b "db' 



a 



dW 

da 



d"¥ 



1 ( dW 



a 



\ * da 



dc 



c 



t « » • 



dc 



• « » I OtS I 
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If we add together the first two of these, and avail ourselves of the property that 
W is homogeneous of degree — 1 , we easily prove that the stationary value of E is 



Jjj «—. -~ r\Jjji 



^ + a 



da 



Since the potential of the ellipsoid must satisfy Poisson's equation 



ixISO 



<PS? clV dW 

ada % bdb cdc 



dW , 7 09 , 09 

a t~ + b — + c 



da 



do 



dc 



__2 
aba 

"9. 



By means of these and two out of the three equations (34). we may eliminate the 
differentials of *P, and writing p for the density find 



or 



111 

"9abc ( — + "xT + 



a 2 ! W l (?- 



6 



\TTp 



(6 s + c 3 ) 



1 i L i 1 



. (35). 



a 



W 



2 



6 



I do not happen to have seen this form for the angular velocity of J acobi's ellipsoid 
in any book. 

It is easy also to show that the stationary value of E may be written 



E = -&M s 



1 1 1 

(6 2 + e % ) ( -jx + To. + 



a? b 2 c 2 



¥ — 2 



J 3 + 6 ,S 

abe 



f 1 



(&» +c a) _+^ + ^)_ 



a 



Z> 3 



6 



c* 



We may now express the potential, say V, of the system entirely in terms of 9 



and a - _— ior 
da 



v 



W 






T a da T S 2 1 da 3i¥ ' ^ 



*3 / 



<W 2« 3 c 3X 



o \ ^ 7 



3i¥ 



+ £*>* (y 2 + z 2 ) , 



w 



d^ fx 2 



r 






* + a l£fe + Vi + * 



We thus verify that V is constant over the surface of the ellipsoid. 

Let g denote the value of gravity at the surface. Then if dn be an element of the 



outward normal, g 



clV 
dn 



Since 



dx 
dn 



px 
a 2 > 



dy py d-z 

dn b 2 ' dn 



V 



,2 ' 



O JL o 
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o> + ?/ + 







/>•* 



£ 



fifa 



6fc ^ 



*2\ 



its / f/ ,</ 

- 4_ ^ - 4. '.'- 
a* T J* T c* 



Now change the notation and write 



3 M dW 



.JL XlvJXJ, 



Now 






h % iv^ — - ] , . h % = F(^ a — 1) 3 c 3 = &V 2 , 



\ 



^ zr: P {y % — £> Q 3 ) 



a 



**" 



da 



9, 



fC j 



,00 



dp 



vo 



(v 



2 _ L±£ 
1 ~/3 



)*(*» - 1) 



/ C) 



M 



K 



1 + /3 W 



oo 



v 






dp 
i + J3Y§- / i 



i 



- i)* " 



i (,) = i, Pi 1 (?) = {s —rzff 



- > 



and 



1?/ fa>) €/ ("<>) - [f .■' ("<>)? ] ^ ^,p {^z^hf(f _ 1)1 ; 



SO UJQciu 



and 



9 



¥ = r $o (*o) <®o ("o) > 



/fc 



cW __ 2 



Pi 1 (v ) Qi 1 (v ) 



* 3 



3i¥" 
^ = "^ P i X ("o) Qi 1 ( v o) 



* « # , lOO|» 



We may note in passing that the condition for a Jacobian ellipsoid (the last 
equation of (34) ) is reducible to the form 



"T~T" ~~T~^ d\h ~ K COt^ y 



> tail 3 i|r 



JQ 



cZi/;. 



On examining the series of functions given in (25), we see that it may be written 

- W (*o) ©a 1 (*o) = Pi 1 (*o) Qi 1 Wo • 

This agrees with M. PoiNCABi's equation (l) on p. 341 of his memoir. 

We will now suppose that the body, instead of being an ellipsoid, is an ellipsoidal 
harmonib deformation of an ellipsoid, which is itself a figure of equilibrium for 
rotation m. 

The addition to E will consist of three parts ; first that due to the mutual 
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energy of the layer of deformation ; secondly that due to the ellipsoid and the 
layer ; thirdly that due to the change in the moment of inertia. 

If a subscript I denotes integration throughout the space occupied by the layer, 
U the potential of the ellipsoid, and dv an element of volume, 



8E^i 



— j- — - + Updv + ^o> 2 (y 2 + s 2 ) pdv. 



If £ denotes the thickness of the layer standing on the element dor, the first 



1 <& 



f^do-jdac 



D 



of these terms is \p 

The value of U + ^ {if + ^) throughout the layer is equal to F — gtf, where 
V is the constant value of U + i^ 2 {}/ + ^ 3 ) over the surface of the ellipsoid, and 
£' is the distance measured along the normal to the element dgdcr of volume. 



Hence 



Updv+ ^co 2 (if + z l ) pdv = p ( V — #£') d£'do- . 

J i J Jo 



Since F is constant and the total mass of the layer is zero, this is equal to 



It follows that 



8E = \p 



1^2 






D 



12 



j g^dcr . 



The axes of the ellipsoid have been chosen so as to make our original E 
stationary, and the further condition to be satisfied is that SE shall be stationary. 
Let us suppose that 

C = pe %> (ji) €f (<!>), 

which expression shall he deemed to include any one of the other types of harmonic. 

Then it is shown in (51) of "Harmonics" that the potential of this layer at the 
surface of the ellipsoid is 



?>M 



I 



7 e %> (v ) <&< (v Q ) %> (/,) €f (4). 



Since the mass of an element is pep l^f (p) (&f (<f>) da-, we have 



1 9 

2H 



m^i = | %> e , m ^ m (vo) 



12 



LW (/*) W WJpdar, 



With the value of g found in (36) 



2P 



giUo- = f f e* P x i (v ) Q^ (v ) 



fflf (p) <&f (4>)J pd<r. 
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Hence 



IV 



FV M Qi 1 ("«) 



W (/0 C' (*)]>?<r. 



In order that the new figure may be one of equilibrium, this expression must oe 
stationary for variations of e. It follows that we must either have e = 0, which 
leads back to Jacobi's ellipsoid, or else 

Pi 1 (* ) Qi 1 (*o) 

This last condition is what M. Poincake calls the vanishing of a coefficient of 
stability. # It shows that if v and /3 satisfy not only the condition for the Jacobian 
ellipsoid, namely, p s T (v ) © 3 T (v ) = P x l (i> ) Q L X (v ), but also this equation, we have 
arrived at a figure which belongs at the same time to two series, and there is a 
bifurcation at this point. The form of the figure is found by attributing to e any 
arbitrary but small value. 



§ 6. The Properties of the Successive Coefficients of Stability. 

Corresponding to each harmonic deformation of the ellipsoid, there is a coefficient 
of stability of one of the two forms 

Pi 1 M Q?~fc) p, 1 W Q? W ' 

These coefficients may be written It/ or K/ according to an easily intelligible 
notation. The Jacobian ellipsoid is defined by v Q , and the question arises as to the 
possibility of the vanishing of the several Jfc's as v Q gradually diminishes from infinity, 
that is to say, as the ellipsoid lengthens. 

An harmonic of the first order merely denotes a shift of the centre of inertia 
along one of the three axes; one of the second order denotes a change of ellipticity 
of the ellipsoid. Since we must keep the centre of inertia at the origin, and since 
the ellipticity is determined by the consideration that the ellipsoid is a Jacobian, 
these harmonics need not be considered, and we may begin with those of the 
third order. 

I shall not attempt to follow M. Poincare in his masterly discussion of the pro- 
perties of the coefficients of stability, f but will merely restate in. my own notation 
the principal conclusions at which he has arrived. 

* i Acta Math./ vol. 7, 1885, p. 321, The factors -J and lj2n + 1 (or l/2i + 1, if i is the degree of 
the harmonic) which occur in his form of the condition are included in my functions. 

t Sections 10 and 12 of his memoir. I have to thank him for saving me from making a serious mistake 
in this portion of my work, 
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1st. The equation 

Pi 1 (") Qi 1 M - W (*0 ©••' (v) or P/ (r) Q/ (v) = 0, (t > 2) 

is not satisfied by any value of v between 1 and infinity, if $ ;' or P/ is divisible by 

v^—*- -r ). It appears from the forms of the functions as given in § 4 of 

" Harmonics " that the P functions are so divisible. These functions appertain to 
the types EES, OOO, OES, EOO, and therefore the ellipsoid cannot bifurcate into 
deformations of these types. 

2nd. The equation has no solution if p* is divisible by (v* — 1)*. We again see 
from § 4 of "Harmonics" that |J/ is so divisible if it is of the types OOS, EOS. 
Hence the ellipsoid cannot bifurcate into these types. The only types remaining 
are EEC, OEO. 

3rd. The equation has no solution if any of the roots of $* (v) = lie outside the 
limits + 1 to — 1 . The only JJ/ of the types EEC, OEC which has all its roots inside 
the limits + 1 to — 1 is the zonal harmonic for which s = 0. 

Hence the ellipsoid can only bifurcate into a zonal harmonic. 

4th. The equation 

Pi 1 Qi 1 -»<©* = (t>2) 

must have a solution between 1 and infinity for all values of i 

It follows from these four propositions that the Jacobian ellipsoid is stable for all 
deformations except the zonal ones, and that as it lengthens it must at successive 
stages bifurcate into each and all the zonal deformations. 

5th. As the ellipsoid lengthens, the first coefficient of stability to vanish is that 
of the third zonal harmonic. This stage is the end. of the stability of the Jacobian 
ellipsoids, and there is almost certainly exchange of stability with the pear-shaped 
figure defined by this harmonic. 

6th. It has not been rigorously proved that there is only one solution of the 
equation 3fc = even in the case where i = 3, but M. Poincaee believes that this 
is almost certainly the case. 

7th. The functions 

W M 1 W (y) ] W (") 1 W (vo) 

or > X or > — or > X or 

P/KU P/WJ P/(v)J P/(v )J 

have always the same sign as v increases from v to infinity, provided that s and t 
are both greater than zero, and i greater than 2. 

The seventh of the preceding propositions renders it easy to determine the relative 
magnitudes of all the It's belonging to a single degree i. 

In what follows I may take the symbols p, © as including also P, Q. 
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ow 



s> — - <C aKi" *^S. 



*.* 



%' ("o) ©i' (*o) - W W ©/ ("o) < - > . 
If we express the ®'s in terms of integrals this becomes 







z/ < = > 



The seventh proposition shows that when s and / are greater than zero, and i is 
greater than 2, all the elements of the integral have the same sign, Hence the 
question is whether. 



»'("o) 

f i s (?) 



< - > 



W K) 



WW ' 



Therefore we have to arrange all the ™ y *! hi descending order of magnitude, 

and shall thereby obtain the non-zonal It's in ascending order, 

T wish first to show that these coefficients may to 'a great extent be sorted by 
considering the inequality. 



§^< = >|^ (,= 1,2,3 



% ] l — — t, ^, o . , ., ?-K 



Suppose, if possible, that whereas, for the ellipsoids defined by /3, i>, v , 

W("o) ^W("o) t Pi'W^ IV (v ) 



WW W(") 



IV <» P'O) ' 



Then there must be some value of /3 for which 

%'{»o)W {») = &• WW fa) 

for all values of v greater than v . 

It is almost obvious that there is no one value of /3 which renders this equation 
possible ; but consider for example the case of s = 2, t = 0. 

Now 

fa 3 W - - Ws (*) + P ? , 2 H f 3 (*) = P s (v) + j89 2 P 8 » (v ). 

Tt we substitute this in the equation we find 

P 3 3 (*o) P 3 W = P 8 9 W P 3 (v ). 

This can only be satisfied by i> = j> , and hence the hypothesis is negatived. 
Similarly the assumption of other values of s and t leads to an impossibility. 
Thus we may consider the P functions in place of the ^ functions. 
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Consider the inequality 

Pi' (Pq) . _ . 1V+ 1 (go) f _ _ 1 „ • 

If the inequality is determined for any value of p, it is determined for all values. 
Now when v is very large 



P^)=^#zr^ *V +1 (") 2 ' ! 



V*. 



2^H- S r 5 ** v 2*ili — s— 1 ! 

Hence our inequality becomes 

(i - *) P/ (v ) > = <P i « +1 (v ). 

This inequality is of the same kind for all values of p . Now P/ (i/ ) involves the 
factor (V 3 — if 8 and P/ + l (v ) involves {v$ — l) i(5 + 1) . Putting therefore p a = 1 + e, 
the left-hand side involves c s * and the right £*<* + 1) . It follows that unless 6' is equal 

O J. 

to { the left-hand side is greater than the right ; but s is necessarily equal to i — - 1 

at greatest. 

Therefore ft' M ^ 1Y+ 1 M 

¥i'(v) Pi'+'O'j' 

Hence K's with smaller 5 are less than those with greater s. 

It remains to discriminate between the two sorts of P-functions which occur in 
ellipsoidal harmonic analysis ; that is to say we must determine 

> = < 



W 00 Pi s « * 



*' 3 



Since the /3 of " Harmonics " is equal to - — — in the present notation, when /3 
and k' are small we have by the formulae of that paper 

W (v) = P> (y) + iK% + J>r 2 (v) + ^^.-zP,-- 2 (v) + . . • , 






p .-' (") + s -~ "%+zPr 2 (*) + V 2 k'V 2 iv- 2 (") + •■ 



—j 



2s ** + <-' v/ ' 2 



When z> is very great and t< f very small ^/ = P/, so it suffices to determine the 
inequality 

%' (v ) > = < P/ ("o) ; 

and this may be considered for any value of v greater than unity. By taking v very 
large and k very small the inequality becomes 

/ l\i 



or .1 j> = <* k . 



i') rn O 
<u 1. zj 
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But k < 1, hence the first sign holds true and 

whence , %f < K,*. 

Thus it follows that for order i 

%" < K l < %r < K, s . . . . <%*< Ki\ 

The order of magnitude of these coefficients is therefore completely determined. 

As confirmatory of the correctness of this result it .may be mentioned that 1 .find 
that when y = 69° 50' and k = sin 7o° 5G' 3 

% { = -1765, Ka 1 = '2990, H/ = '4467, K 3 3 = 9 4550 3 »/ = '5719, K 3 3 = '5876. 

When y == 75° and k = sin 81° 4 ; (another Jacohian ellipsoid) the numbers run 
•130, '224, -460, '465, '604, '614. 

We see that for the harmonics of higher order the ellipsoid is more stable than it 
was and for those of lower order less stable. 



§ 7. The critical Jacohian Ellipsoid. 

From a number of preliminary calculations I saw reason to believe that the critical 
ellipsoid would be found within the region comprised between y = 69° 48' and 
69° 50', and sin™ 1 k == 73° 52' and 73° 56'. 

If we write 

/./ • ™i \ E / . /e 1 ' sin 3 7 cos 3 7V /oTn 1?x « 3 sin 7 cos 7 (1 + ic 3 sin 3 7' 
/ (y, sm i K) = 7 Jl+ T~T^7sT; "" ( 21 "" L ) 



*; " \ 1 — k" sm v 7 / v y /e ,J (1 -— a:" sur 7)* 



■where the amplitudes of E and F are y and their moduli /c, the existence of the 
Jacobian ellipsoid is determined by 



f(y, sin l k) — 0.* 



The coefficient of stability is 



*fr / en\,-l \ — 1 ftW^aW 

^( y ,sm ^-I-pT^-qT^)' 

The formulae for computing lt 3 are given in § 4. 
The values of E and F are from Legendre's tables . 

* See <lioy. Soc. Proc., 5 vol 41, p. 323^ where the formula, m reduced to a form convenient for 
computation. 
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Now I find 

/(69° 48', 73° 52') = + '000191 ; /(69° 50', 73° 52') = + '001319. 
/(69° 48^ 73° 56') = - '001186 ; /(69° 50', 73° 56') = - '000031. 

» 3 (69° 48\ 73° 52') = + '001058 ; %(69° 50', 73° 52') = - '000885. 
% (69° 48', 73° 56') = + '000655 ; lt 3 (69° 50', 73° 56') = - '000765. 

By interpolation we get the following results : — 
The Jacobian ellipsoid is given by 

( y — 69° 48') — '59642 (sin™ 1 k ~~ 73° 52') + '33091 = 0. 

The vanishing of the coefficient of stability is given by 

(y — 69° 48') + '041625 (sin"" 1 k - 73° 52') - 1/0890 - 0. 

In these equations the minute of arc is the unit. 
Solving them I find 

y = 69° 48 / *997 = 62° 49''0, 
sijr 1 * = 73° 54 ; '225 = 73° 54''2. 
With these values I find that the three axes a, &, c, where ate = a :J are 

'- = '650659, 

a 

" = -814975, 

a 

-- = 1*885827. 

The last place of decimals in these is certainly doubtful. 
The formula for o> 3 is given in (35). 

Now V = \ f (v ) ©o (*o), * = <* sin y, % (v ) ©, (y„) = kY. 

lb 

Then since a = c cos y, 6 = cA, 

2ttp ~~~ 1 4- A 3 — - — 

In this formula, F, y, A must correspond with values interpolated, amongst those 
used in obtaining the solution. 
From this I find 

?".. = -1419990 = -14200. 

'2lTp 
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In the paper on the Jacobiaii ellipsoid referred to above the moment of momentum 
is tabulated by means of ji, where the moment of momentum is (f ?rp)ia 5 /x, The 
formula for [m is given in (25) of that paper, and, modified to suit the present 
notation, is 



/x = ^(Acosy)-*(l + A 2 ) 



w V?; 



5 



4 



irp 



For the critical ellipsoid I find ft = '389570. 

The following table gives the numerical values for a number of Jacobiaii ellipsoids, 
beginning with the initial one and terminating just beyond instability. The last line 
;ives the corresponding values for the critical ellipsoid. 



8 



Jacobi's Ellipsoids. # 
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57 . . 


* 9 * 


34f 
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60 . . 


* C * 


49 


^7 


40 54 


ft r* 

bo . . 


* » * 


64 


19 


54 46 


70 . . 


W • • 


74 


12 


64 43 


69 49 


« • 


73 


04: 


64 24 



a /a. 



•6977 

•697 

•696 

•6916 

•6765 

•6494 

65066 



/>/a. 



1-1972 ! 

1-179 | 

1*0454 I 

•9235 ! 

•8111 i 

F 

•81498 I 



eya, 



1-1972 

1-216 

1-279 

1-3831 

1-6007 

1-899 



1-88583 



W 3 /27T/J, 



'18712 

18706 

•186 

•1812 

•1659 

1409 

14200 



^x» 



^30375 

•304 

•306 



•3134 
•3407 

•3920 



•38957 



* I have been criticised with respect to my paper on Jacobi's ellipsoid, from which these results are 
extracted, by M. S. Kruger (Nieuw Archief voor Wiskunde, Tweecle Keeks, Derde Deel and 
' Ellipsoidale Evenwic fits vor men/ &c, Thesis for Degree of Doctor, Leiden, J, "YV. van Lceuwen, 
Hoogewoerd 89, 1896), because I wrote it in ignorance of certain previous work, especially of a paper 
by Plana (' Ast. Nadir./ 36, n. 851, c. 169), But I cannot but congratulate myself on my ignorance, 
since it appears that Plana gave a number of numerical results which were wholly wrong. A knowledge 
of that paper would no doubt have caused me much further trouble. 

My paper gives a number of solutions of the problem which I believe to be correct. Unfortunately 
the methods of the paper are clumsy, and there are several mistakes. The formula for oj 2 used in this 
present paper, is much better than that given there. 

The complicated formula on p. 325 is susceptible of reduction to a simple form, for on substituting for 
y its approximate form (i) we have simply 

y - S = |-k 2 sin 8 cos 8, 

where 8 - 54° 21 ; 27". 

The hnal numerical result was, however, nearly right, for I now find 

sin a a = 10^260821 sin (y - 8), 

whereas I had '9266528. The sin a is the same as the k used here, 

The formula at the top of p. 326 which is reproduced as (22) on p. 828 is, I think, illusory, for if in the 
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In order to determine the question as to whether or not it is possible that % % = 
should have another solution than that found in the next section, I have computed 
the value of this coefficient for the Jacobian ellipsoid y = 75°, k = sin 81° 4'*4, and 
find it to be — 6*627. From the manner in which the numbers in the computation 
present themselves, it is obvious that for more elongated ellipsoids Jfc 3 will always 
remain negative, and will become numerically greater, I have therefore not thought 
it necessary to seek for an algebraic proof that there is no second root of the 
equation. 

Very long Jacobian ellipsoids tend to become figures of revolution, and the coeffi- 
cients of stability tend to assume the forms 



*V (") Qi 1 (v) ' 



The forms of these functions are well known, and I think that fair approximations 
to the incidences of the successive figures of bifurcation might be derived from the 
vanishing of this expression. 

For example 

Pi 1 W Qi 1 (*) = * [* - V* - i) % ^l 

P 4 (*) Q*(") = 6 J 4 (35^ - 31V + 3) log (^1J - iv (21*» - 1 1 j (35i>* - 30*» + 3) 

I have not, however, attempted to solve the equation found by equating these two 
expressions to one another. 

Even when i = 3 and y = 69° 49' (the critical Jacobian) this rough approximation 
makes the coefficient of stability very small, but it is to be admitted that P^ Qj 1 and 
P 3 Q 3 differ very sensibly from P^ (v) Q^ (v) and |3 8 (v) <B 3 (v), although in such a 
way that the errors compensate one another. 

first term we pub y = S 4- \k~ sin 8 cos $ (as is clearly allowable in approximation) the term with coeffi- 
cient k~ or sin 2 a disappears. This shows that it was necessary to proceed in the approximation as far as 
k 4 in order to obtain a result. 

The methods of approximation adopted on pp. 326-7 are of doubtful propriety, but will, I think, lead to 
nearly correct results. There is, however, a mistake towards the bottom of p. 327 which runs on to the 
end. M. Kruger correctly points out that the second line of formula (24) p. 329 should run 



J cos 2 a 



-i- log, cot (£tt - Jy) . ( J # + 3 tan 2 y + tan 4 r ) - - 1 / ~ ¥ tan 2 7 ~ t tan 4 7 
.sin 7 



Lastly, on p. 335, line 13, for = 0*3573, read C = 0-5379; and on p. 336, line 7, for 1*3573, read 



1-5379; and for- - 1-696, read-- 4-65. 

a a 
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§ 8. The pear-shaped Figure, of Equilibrium. 

By (21) the normal displacement hn for the third zonal harmonic deformation may 
he written 

« z [#' 2 £ 2 — qhP — (3 — 4-q^) ?/ 3 — c 3 $ V 3 sin 9 7] 

on = £ -~ 



r y»(l _ ^ R i n s 7 ) (;^/ C os 4 7 + ? /VA 4 _ + S 2 )* ? 



subject to the condition 



9 o 

— -J- ^ = c~, 



The expression has been arranged so that when x = ?/ = 0, 2 = 6\ we have Sn = e. 
Hence + e and .— 6 are the normal displacements at the stalk and blunt end of the 
pear respectively. 

In the section y = 0, this may be written 



q ~ c"(rr — #- snv 7)-' 

The nodal points are given by ~ = ± 9 = ± '758056, 

In the section # = 0, since k j = <f~ — ^ > it may be written 



Sn 



A (4 — 5^ 3 ) # (#% 3 — c 3 ^ 3 ) 



q ,J etc* (a j — /CZ" Bur 7)- 



At jy 

The nodal points are given by : = + ~ = + '788986. 



1 

c ~ tc 

The section % = is obviously another nodal line for all sections. 

By means of these formulae it is easy to compute the normal displacements from the 
surface of the critical Jacobian. 

The figure opposite showing the three sections ,^ = 0^ y = 0, 2 = 0, is drawn 
from these formulae, the dotted line being the critical Jacobian and the firm line the 
pear. The scale of the normal displacements is, of coimse, arbitrary. 

Comparison with M. Poinoare's sketch shows that the figure is considerably longer 
than he supposed. 

In this first approximation the positions of the nodal lines are independent of the 
magnitude of e 9 and they lie so near the ends that it is impossible to construct an 
exaggerated figure, for if we do so the blunt end acquires a dimple, which is absurd. 
It might have been hoped that such an exaggeration would afford us some idea of the 
mode of development of the pear. 
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M. Schwarzschild has remarked* that it is not absolutely certain that the 
principle of exchange of stability holds with reference to this figure, and that we 
cannot feel absolutely certain that the pear is stable unless we can prove that the 
moment of momentum is greater than in the critical Jacobian. 

With reference to this objection, M. Poincare writes to me as follows : — 

" Faisons croitre le moment de rotation, que j'appellerai M. Deux hypotheses sont 
possibles. 

" Ou bien pour M < M (the moment of momentum of the Jacobian), nous aurons 
une seule figure, stable, a savoir l'ellipsoide de Jacobi, et pour M > M trois figures, 
une instable, Tellipsoide, et deux stables (d'ailleurs egales entre elles), les deux figures 
pyriformes. 

" Ou bien pour M < M , nous aurons trois figures d'equilibre, deux pyriformes 
instables, une stable, l'ellipsoide, et pour M < M une seule figure instable, l'ellipsoide — 
auquel cas la masse fluide devrait se dissoudre par un cataclysme subit. 

" II y a done a verifier si pour les figures pyriformes, M > ou < M ." 

It seems very improbable that the latter can be the case ; but this opinion is not a 
proof. 

Since o) 2 is stationary for the initial pear, a small change in the angular velocity 
will certainly produce a great change in the figure of the pear. If this investigation 
has, in fact, its counterpart in the genesis of satellites and planets, it seems clear that 
the birth of a new body, although not cataclysmal, is rapid. 



§ 9. Summary. 

It is possible by the methods explained in my previous paper on " Harmonics" to 
form rigorous expressions for the ellipsoidal harmonics of the third degree. Accordingly 
in § 1 I proceed to form those functions. In § 2 the notation is changed with a view 
to convenience in subsequent work, and for the sake of completeness the harmonics 
of the first and second degrees are also given. In § 3 the corresponding solid 
harmonics are expressed in rectangular co-ordinates x 5 y 9 z. In § 4 I find the 
Q-functions, the harmonic functions of the second kind, and express the results in 
terms of the elliptic integrals E and F. It appears that both the P- and Q-functions 
of the third degree of harmonics occur in three pairs which have the same algebraic 
forms, and that in each pair one of them only differs from the other in the value of a 
certain parameter. There is, lastly, a seventh function which stands by itself; this 
last corresponds to the solid harmonic xyz. 

In § 5 the equations for Jacobi's ellipsoid are determined by the consideration that 
the energy must be stationary, and the superficial value of gravity is found in terms 
of the appropriate P- and Q-functions. I then proceed to find the additional terms 



* a 
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in the energy when the mass of fluid is subject to an ellipsoidal harmonic deformation. 
This section is a paraphrase of M. Poincare's work, but the notation and manner of 
presentation are somewhat different. The additional terms in the energy are shown 
to involve a certain coefficient, which is called by M. PousrCARE a coefficient of 
stability. It is clear that when any coefficient vanishes we are at a point of bifurca- 
tion, and the particular Jacobian ellipsoid for which it vanishes is also a member of 
another series of figures of equilibrium. 

In § 6 the principal properties of these coefficients, as established by M. Poincare, 
are enumerated. He has shown that the ellipsoid can bifurcate only into figures 
defined by zonal harmonics ; that it must do so for all degrees, and that the first 
bifurcation occurs with the third zonal harmonic. The order of magnitude of the 
coefficients of the several orders and of the same degree is determined. A numerical 
result seems to indicate that as the ellipsoid lengthens, it becomes more stable as 
regards deformations of the third degree and of higher orders, and less stable as 
regards the lower orders of the same degree. 

In § 7 the numerical solution of the vanishing of the coefficient corresponding to 
the third zonal harmonic is found, and it is shown that the critical ellipsoid has its 
three axes proportional to '65066, '81498, 1*88583, and that the square of the angular 

of the axes of the Jacobian ellipsoids from their beginning on to instability at this 
critical stage. The nature of the formula for the third zonal coefficient of stability 
seems to show that it can only vanish once — a point which it appears that 
M. Poincare found himself unable to prove rigorously. 

A suggestion is made for the approximate determination of the bifurcations into 
the successive zonal deformations, but no numerical results are given. 

In § 8 the nature of the pear-shaped figure is determined numerically, and the 

* 

reader may refer to the figure above, where it is delineated. It will be seen to be 
longer than was shown in M. Poincabe's conjectural sketch. 

If, as M. Poincare suggests, the bifurcation into the pear-shaped body leads 
onward stably and continuously to a planet attended by a satellite, the bifurcation into 
the fourth zonal harmonic probably leads unstably to a planet with a satellite on 
each side, that into the fifth to a planet with two satellites on one side and one on the 
other, and so on. 

The pear-shaped bodies are almost certainly stable, but a rigorous and conclusive 
proof is wanting until the angular velocity and moment of momentum corresponding 
to a given pear are determined. To do this further approximation is needed. 



velocity is given by ^^ = "14200. A short table is also given showing the march 
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